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Abstract-An approximation method using to estimate the 
influence of the uncertainties of the neutron flux characteristic 
parameters in the irradiation positions on the NAA results using 
k -standardization technique was presented. These parameters 
are the epithermal reactor neutron spectrum shape-factor a, the 

effective resonace energy E r for a given nuclide and the thermal 
to epithermal neutron flux ratio f. The method is applied to 
estimate the effect of the uncertainties in the determination of a, 

E r and f on final NAA results for some irradiation channels of 

the Dalat reactor. The obtained results shows that presented 
method is suitable in practical use for the estimation of the 
errors due to the uncertainty of the neutron flux characteristic 
parameters at the irradiation position in reactor. 

Keyword-ko-standardization Technique, Error Propagation 
Function, Neutron flux Characteristics, Dalat reactor 

I. INTRODUCTION 

Since the k -standardization method was introduced in 
NAA [1], it has been broadly applied in the reactor in the 
world. The fundamental concept of k -method was being 
elaborated previously in great detail [1-3]. The concentration 
of a element in the k -method is calculated by: 



p(ppm) 



N p/t m 

SDCW 1 f +Qo («) e p 



(1) 



N p/t m 

SDCw 



k f + Qo(«) e p 



with k in Eq.(l) defined as: 

M*9c y 
° M6*a V 
InEqs.(l)and (2): 

M - atomic mass; 

9 - isotopic abundance; 

Go - 2200 ms" 1 (n,y) cross-section; 

y - absolute gamma-intensity; 

N p - peak area corrected for pulse losses; 

W - sample weight in gram; 

w* - comparator weight in microgram; 



(2) 



S = l-exp(-Xt ilr ); t^ - irradiation time; X - decay 

constant; 

D = exp(-Xt d ); t d - decay time; 

C = [l-exp(-^t m )]/Xt m ; t m - measuring time; 

f - thermal to epithermal neutron flux ratio; 

Qo(a) = I (a)/a ; Io(cx) - resonance integral 
corrected for a non-ideal epithermal neutron flux distribution 
( assumed l/E l+a ); 

8 P - detector's efficiency; 

When the epithermal neutron flux distribution deviates 
from ideality, i.e. it does not follow the 1/E-law, Qo(oc) of 
nuclide i can be writen by: 

Qoi(a)=(Q 0i - 0.429)/( E~ ) a +0.429/[(2a+l)(0.55) a ] 

(3) 

with a - neutron spectrum shape factor deviating 
from the 1/E-law, independent of neutron energy and \a\ «1; 

Eri - effective resonace energy of nuclide i; 

The asterisks in Eqs.(l) and (2) refers to the comparator, 
which is suitable for coirradiation with the sample; in most 
case, Au is used as a comparator. The k -factors to Au for 
interested isotopes in NAA were experimentally determined 
and tabulated in [4] with an accuracy better than 2% (average 

~ 1%). The relevant nuclear data as Q oi and Eri can be found 
in a tabulated form or in a computer library, a, f and e p must 
be determined by experiment and they depend on specific 
irradiation channel and detector, which are used in practice. 
The detector's efficiency (e p ) can be determined with 
anuncertainty about 2%; but the uncertainty of a can be more 
than 10%, even bigger, depend on the irradiation channels in 
reactor. Since the term [f + Q (a)]/[f + Qo(a)] in Eq. 

(1), it is clear that an additional parameter, Eri , should be 

considered, because the uncetainties of Eri of some nuclides 
are about 20% [4,5]. 

The accuracy and the applicability of the k - 
standardization method were detailly presented in [5] by F. De 
CORTE et. al.. In [6], J.OP De BEEK evaluated the effect of 

errors of a and Erion the results in terms of concentration, 
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based on the Au comparator; in that Qoi(oc) was 
approximated by : 

Qoi(a)«Q 0i (Eri)- a (4) 

However, with this approximation, it led that some results 
in [6] have to be put to disscusion (see below). 



In this work, we carry out an approximation method to 

evaluate the effect of errors of a and E r i on the NAA results 
in the k - standardization method. The obtained results 
showed that the approximate method in this work is 
acceptable with confident accuracy. 



TABLE I THE VALUES OF &i FOR THE INTERESTING NUCLIDES IN NAA 



Target 


Formed 


a. 




Target 


Formed 


ill 


3j 


Nuclide 


nuclide 


with a<0 


with a>0 


Nuclide 


nuclide 


with oc<0 


with oc>0 


23„T 

Na 


24* t 

Na 


0,524987 


0,301907 


1 16n 

Sn 


1 17nif, 

Sn 


0,996203 


0,992120 


26» * 

Mg 


27» * 

Mg 


0,600414 


0,260458 


122„ 

Sn 


123m„ 

Sn 


0,958568 


0,907080 


27 * i 

Al 


28 * i 


0,632092 


0,355905 


124n 

Sn 


125mo 

Sn 


0,996581 


0,993233 


CI 


CI 


0,618735 


0,340898 


Sb 


122m Sb 


0,99639 


0,991949 


41 -ir 
K 


42t^ 

K 


0,744902 


0,510677 


123 ni 

Sb 


124m Sb 


0,993951 


0,988375 


Sc 


Sc 


0,225177 


0,13871 


127t 

I 


128 j 


0,991934 


0,984365 


Tl 


51,™- 

Ti 


0,61391 


0,300434 


Cs 


134m,-. 

Cs 


0,993476 


0,980868 


5|., 

V 


52* , 

V 


0,484063 


0,253348 


130™ 

Ba 


131m™ 

Ba 


0,991752 


0,983823 


50 Cr 


51 Cr 


0,444713 


0,235221 


132 Ba 


133mg a 


0,961113 


0,921759 


55» » 

Mn 


56» » 

Mn 


0,768236 


0,591296 


138™ 

Ba 


139™ 

Ba 


0,721393 


0,434897 


58™ 

5 Fe 


59™ 

39 Fe 


0,744577 


0,553371 


1 39 T 

y La 


La 


0,824392 


0,703294 


Co 


60*-. 

Co 


0,899648 


0,792867 


140 .-i 

Ce 


141^ 

Ce 


0,696096 


0,429943 


Ni 


Ni 


0,603571 


0,327784 


1 49 ^ 

Ce 


Ce 


0,798767 


0,600323 


63^i 

Cu 


Cu 


0,786528 


0,594356 


141™ 

Pr 


142nipj. 


0,84513 


0,710838 


Cu 


Cu 


0,768558 


0,578803 


146* T i 

Nd 


147*. t i 

Nd 


0,884097 


0,750436 


Zn 


65 rj 

Zn 


0,879728 


0,716396 


148* T i 

Nd 


149-KT 1 

Nd 


0,956282 


0,908206 


Zn 


69m Zn 


0,928832 


0,842012 


1 50tl.T 1 


151».T 1 

'Nd 


0,982296 


0,962764 


71 Ga 


Ga 


0,967117 


0,932900 


1 S9.™. 

L2 Sm 


i 


0,995457 


0,985053 


75 * 

As 


76 * 

As 


0,984389 


0,968602 


154o 

Sm 


155o 

Sm 


0,949148 


0,899154 


74 r, 

Se 


Se 


0,982325 


0,966392 


Eu 


Eu 


1,002410 


0,972428 


Br 


80m Br 


0,984481 


0,963835 


158/-< i 

Gd 


159/-< j 

Gd 


0,993709 


0,987888 


Br 


82m™ 

Br 


0,988802 


0,976569 


160^. i 

Gd 


161 > 

Gd 


0,941083 


0,869589 


Rb 


86m Rb 


0,985106 


0,962061 


159™ 

Tb 


160mi 

Tb 


0,991765 


0,983428 


Rb 


88™, 

Rb 


0,990543 


0,978106 


164™ 

Dy 


165m Dy 


-0,59612 


-0,13894 


Sr 


Sr 


0,984748 


0,963917 


165 TT 

Ho 


166 t t 

Ho 


0,989663 


0,976313 


Sr 


87m^ 

Sr 


0,945262 


0,87112 


170™ 

Tir 


171™ 

Er 


0,950799 


0,903065 




90mY 


0,963615 


0,891243 


Tm 


170rr, 

Tm 


1,004700 


0,991756 


94 Zr 


95 Zr 


0,957566 


0,87700 


174 Yb 


175 Yb 


0,357880 


0,221067 


% Zr 


97 Zr 


0,999115 


0,997943 


176 Yb 


177 Yb 


0,908488 


0,809688 


93 Nb 


94m Nb 


0,969700 


0,928380 


175 Lu 


176m Lu 


0,996032 


0,991774 


98 Mo 


"Mo 


0,995883 


0,990853 


174 Hf 


175 Hf 


0,759887 


0,609407 


100 Mo 


101 Mo 


0,988295 


0,970848 


179 Hf 


180mj^ 


0,990364 


0,980132 


96 Ru 


97 Ru 


0,991702 


0,978811 


180 Hf 


181 Hf 


0,913130 


0,837096 


102 Ru 


103 Ru 


0,938892 


0,877634 


181 Ta 


182m Ta 


0,997216 


0,992780 


104 Ru 


l05 Ru 


0,982706 


0,958954 


186 W 


187 w 


0,988597 


0,977570 


103 Rh 


104m Rh 


1,29702 


1,11028 


185 Re 


186 Re 


1,014020 


0,998171 


108 Pd 




0,993450 


0,987460 


187 Re 


188m R£ , 


0,958039 


0,921987 


nop d 


lllm pd 


0,989017 


0,971347 


190 Os 


191m Os 


0,891403 


0,800890 


107 Ag 


108 Ag 


0,934276 


0,880183 


192 0s 


193 0s 


0,908020 


0,831151 


109 Ag 


■10m Ag 


1,00012 


0,990589 


m b 


194 Ir 


1,049670 


1,015320 


114 Cd 


115 Cd 


0,994496 


0,988007 


198 Pt 


199mp t 


0,987525 


0,974821 


113 In 


114m In 


0,999628 


0,993641 


197 Au 


198 Au 


1,001300 


0,990335 


115 In 


116m In 


1,07891 


1,03788 


1% Hg 


197m Hg 


0,493779 


-0.32989 


112 Sn 


113m Sn 


0,995628 


0,991087 


238 U 


239 u 


1,0004 


0,99725 



II. BASE OF APPROXIMATION 
In the approximation of J. OP De BEEK, it is good for the 
nuclides having Q oi >1, but is not for the nuclides with Q oi <1. 

Indeed, in case of 45 Sc(n,y) 46 Sc with Q = 0.44, E r = 5130 eV, 
if a = 0.1, the value of Q (a) from Eq. (4) of J. OP De BEEK 
is 0.18, but from Eq. (3) it is 0.383. The difference is a factor 
of 2. If a = 0.2, Q (a) is 0.08 and 0.345 fromEq. (4) and Eq. 
(3), respectively. The difference is more than four times (!). 
Correspondingly, some results in [6] are not accurate and 
needed to be discussed. 



As we know, a value is smaller than unity in absolutte 
value. In practice, in irradiation channels of reactor, absolute 

value of a is less than 0.2 (in most cases, \(%\ < 0.1 and this 
condition is satisfactory in reactor core). For this reason, 
Qoi(a) from Eq. (3) can be written in a simpler appximative 

expression. Due to \ a \ « 1, in [7,8], we suggest substituting 

Q j (a) from Eq. (3) by the following approximated formula: 

Qo I (a)«Qoi(Eri)" a i a or 

Qo, («) * Q 0i exp(-a, (lnE~)«) (5) 



IJNESE Vol.1 No. 1 2011 PP.22-27 www.ijnese.org ©World Academic Publishing 



-23- 



The International Journal of Nuclear Energy Science and Engineering 



IJNESE 



Where aj is constant for each nuclide and determined by 
fitting the values of Q i(a), which are calculated from Eq. (3) 

in range \ a \ < 0.2 , then fitting according to function (5) (see 

[7,8]). Note that, a ; of each nuclide depends on the sign of a. 
The values of a; for the interested nuclides in NAA are given 
in Table 1. Seeing the Eq. (5), it differs to Eq. (4) of J. OP De 
BEEK by a correctional coefficient a;. However, it can be 
used good for all nuclides with uncertainties of the 

calculated Q i (a) less than about 5% for the nuclides having 

Q D1 <1 and less than about 2% for Q oi >1 with |«| < 0.2. 

Indeed, we carried out a survey of the ratios of Qoi(a) 
calculated from Eq.(5) (in this work) and Eq. (4) (of J. OP De 
BEEK) to Eq. (3) (accurate expression) for Q i from 0.44 
( 46 Sc) to 248 ( 97 Zr) with a=-0. 1 and the results are presented 
in Fig. 1 . 



TABLE II THE VALUES OF Q oi (a) CALCULATED FROM THREE 
EXPRESSION Eq.(3), Eq. (4) AND Eq. (5) 
WITH a IN INTERVAL [-0.2,0.2] FOR 45 Sc(n ,y) 46 SC 



Value 


Qoi(a) 


Qoi(a) 


Qoi(a) 


Value 


Qoi(a) 


Qoi (a) 


Qoi (a) 


of a 


from 


from 


from 


of a 


from 


from 


from 




Eq.3 


Eq.5 


Eq.4 




Eq.3 


Eq.5 


Eq.4 


-0,20 


0,7072 


0,6850 


2,429 


0,02 


0,4263 


0,4287 


0,371 


-0,18 


0,6629 


0,6554 


2,048 


0,04 


0,4139 


0,4178 


0,313 


-0,16 


0,6242 


0,6270 


1,726 


0,06 


0,4026 


0,4071 


0,264 


-0,14 


0,5905 


0,5998 


1,455 


0,08 


0,3923 


0,3967 


0,222 


-0,12 


0,5607 


0,5738 


1,227 


0,10 


0,3828 


0,3865 


0,187 


-0,10 


0,5345 


0,5490 


1,034 


0,12 


0,3741 


0,3767 


0,158 


-0,08 


0,5112 


0,5252 


0,872 


0,14 


0,3661 


0,3671 


0,133 


-0,06 


0,4905 


0,5025 


0,7346 


0,16 


0,3587 


0,3577 


0,112 


-0,04 


0,4718 


0,4807 


0,619 


0,18 


0,3518 


0,3486 


0,095 


-0,02 


0,4551 


0,4599 


0,522 


0,20 


0,3455 


0,3396 


0,080 



Notice: Eq.3: true expression Eq.5: Expression in this work Eq.4: Expression 
in [6] of J. Op De Beek 



2.5 1 




TABLE m CHARACTERISTICS OF ISOTOPES USED IN THE 
CALCULATION IN THIS WORK 



Nuclide 


Q =I,)/CFo 


E r (eV) 


45 Sc 


0,44 


5130 


59 Co 


1,993 


136 


94 Zr 


5,05 


6260 


1S6 W 


13,7 


20,5 


197 Au 


15,7 


5,65 


9s Mo 


53,1 


241 


96 Zr 


248 


338 



III. RESULTS AND DISCUSSION 

The absolute uncertainty in p can be calculated from the 
uncertainties of the variables (denoted Xj) which determine p 
in Eq. (1): 



t Ratio of Q i(a) calculated by Eq.(5) in this work to Eq.(3) 

Ratio of Q oi (a) calculated by Eq.{4) of I. OP De BEEK to Eq.(3) 

Fig 1. Survey of the ratios of Qoi(oc) calculated from Eq.(5) (in this work) and 
Eq. (4) (of J. OP De BEEK) to Eq. (3) (accurate expression) for different Qoi 
with a = -0. 1 

Clearly, the approximated expression in this work is better 
than one of J. OP De BEEK. Moreover, the calculated Qoi(a) 
from three expression Eq.(3), Eq. (4) and Eq. (5) 
for 45 Sc(n,y) 46 Sc presented in Table 2 and the another nuclides 
presented in papers [7,8] also confirm the above conclusion. 

From Table 1, it shows that coefficients a; of nuclides 
having Q i >1 are close to unity, but a; of the nuclides having 
Qoi <1 differs more than unity. Therefore, the approximation 
of Eq. (4) in paper [6] is only acceptable for the nuclides 
having Q oi >1, but for the nuclides having Q j <1, it is not 
reliable. 



E 



dp 

v dx jy 



(6) 



where da/dxj are the corresponding partial derivatives. 

TABLE IV CHARACTERISTICS OF IRRADIATION CHANNELS 



Channel 


a 


f 


Channel 7-1 (Dalat reactor) 
Neutron trap (Dalat reactor) 
Channel 17 (Thetis reactor)[9] 


-0.044 + 0.004 
-0.03 1± 0.004 
-0.028 


14.2 + 0.5 
33.0±0.5 
15.0 



According to the customary error propagation theory, the 
error propagation functions can be written as: 



Z p ( Xj ) 
and relative error is: 



clp 
P 



v x i J 



dp *j 
5x p 



(7) 



In this work, we use the approximation expression; Eq. (5), 

to evaluate influence of the uncertainties of a, f and E ri On 
the final element concentration in k -method in the channels; 
7-1, neutron trap of Dalat reactor (Vietnam) and channel 17 of 
THETIS reactor (Belgium) for the 

nuclides; 45 Sc, 59 Co, 94 Zr, 186 W, 197 Au, 98 Mo, 96 Zr. We choose 
these nuclides, because they differ considerably in Q oi and 

Eri values. The numerical data of concerning isotopes and 
irradiation channels used in this work are summarized in 
Tables 3 and 4. 



s D (x i ) = Z D (x i )^ 



(8) 



1 ) Influence of Uncertainty of E ri on NAA Results 

From Eq. (8), the uncertainty of the concentration (p) in 
k -method due to the uncertainties of the effective resonace 
energies can be written by: 
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Ap 
P 



\ AE ri 

' En 



(9) 



Using Eq.(7) for the effective resonance energy of the 
nuclide i, we obtain: 



a 



a iQoi 



V 



Q 0i + 



E ri )' 



(10) 



The values of calculated Z p (E r i) for choosen nuclides 

are presented in Table 5. The effect of the effective resonance 
energy on NAA result include the uncertainties of the 
effective resonance energies of analytical and comparator 

nuclides. In this case, Au used as comparator with Ea« of 
5,65 eV and uncertainty of 7.1% from [5], the contribution of 

the uncertainty of Ea« to the error of NAA result in 
channels 7-1, neutron trap of Dalat reactor and channel 17 of 
THETIS reactor is 0,17%, 0,077% and 0,13%, respectively. 
Clearly, the effect of the uncertainty of the effective resonance 
energy of Au is negligible and can be overlooked in the 
evaluation. 



TABLE V CALCULATED RESULTS OF 
NUCLIDES 



Z p( E *) 



FOR CHOOSEN 



Nuclide 


7-1 channel 
(Dalat reactor) 


Neutron trap 
(Dalat reactor) 


Channel 17 

(THETIS 
reactor) 


45 Sc 


0,000316 


0,000102 


0,000185 


59 Co 


0,005765 


0,001939 


0,00329 


94 Zr 


0,01430 


0,005304 


0,008004 


185 W 


0,02278 


0,01025 


0,01379 


""Au 


0,02396 


0,01106 


0,01467 


98 Mo 


0,03620 


0,02163 


0,02244 


% Zr 


0,04601 


0,03065 


0,02671 



The analysis for 94 nuclides used in NAA showed that the 
uncertainties of their effective resonace energy are from to 
20%, except 75 As (34%) [4]. In this measure, we are able to 
realize that the effect of them on NAA result is also negligible. 

For example, 45 Sc ( E r =5130 eV, AE r =17%) and 95 Zr 

( Er =338 eV, AE r =2.1%), the contribution of the 
uncertainty of the effective resonance energy to the error of 
NAA result in three above channels is less than 0.01% 
for 45 Scand 0.1%for 95 Zr. 

In epicadmium neutron activation analysis (ENAA), the f- 
term in Eq. (10) should be omitted. The error propagation 

function of E r i can be written: 



Z p( E ") = 



(11) 



The calculated results of 



Z p ( 

197 F 



E ri ) 



for for the 

nuclides; 45 Sc, 59 Co, 94 Zr, 186 W, l97 Au, 98 Mo, 96 Zr in ENAA 
are carried in Table 6. In this case, the error propagation 
function is higher than in the one of irradiation without 
cadmium. Generaly speaking, a; <1 and if a « 1, the 

contribution of AEri to the error of NAA result for almost 



analytical nuclides is less than 1 % and can be omitted in the 
calculation. 

TABLE VI CALCULATED RESULTS OF Zp( E rf) FOR THE NUCLIDES 

IN ENAA 



lViicliHp 


7-1 channel 
(Dalat reactor) 


Npntrnn tran 
(Dalat reactor) 


Channel 17 

(THETIS 
reactor) 


45 Sc 


0,00991 


0,00743 


0,00631 


59 Co 


0,03958 


0,02969 


0,02519 


' 4 Zr 


0,04213 


0,03160 


0,02681 


185 W 


0,04350 


0,03262 


0,02768 


""Au 


0,04406 


0,03304 


0,02803 


98 Mo 


0,04382 


0,03286 


0,02788 


96 Zr 


0,04396 


0,03329 


0,02797 



2 ) Influence of Uncertainty of a on NAA Results 

Also from Eq. (8), the uncertainty of p due to 
uncertainty of a can be written: 



the 



\ Aa 
t) — 



L=Z p (a 
p a 

and error propagation function of a 

f I \ 

a s Qm In E ri 



Z p (a) = 



-a 



a iQoi ln 



Qoi +f E 



+ 



\ a i a 



(12) 



13) 



The values of the error propagation function of a in the 
channels; 7-1 and neutron trap of Dalat reactor (Vietnam) and 
channel 17 of THETIS reactor (Belgium) for the 
nuclides; 45 Sc, 59 Co, 94 Zr, 186 W, 197 Au, 98 Mo, 6 Zr were shown 
in Table 7. From Table 7, for the nuclides having Qo<Qoa u in 
three these channels, the contribution of the uncertainty of a 
to the error of NAA result is not significant, about less than 
1%. But for nuclides having Q »Q au, this effect is 
noticeable. For instance, in channel 7-1 of Dalat reactor (a =- 
0.044, Aa =12% [7,8]), the contribution of the uncertainty of 
a on the error of result of 45 Sc (Q =0.44) is 0.42%, but 
for 99 Mo and 96 Zr is 1,36% and 2.4%, respectively. As a 
comment, for RNAA using 197 Au comparator, the systematic 
effect for a value up to 0.1 is practically negligible for all 
nuclides with a low enough Q value (e.g. 45 Sc, 59 Co, 58 Fe, 
ect.). On the other hand, for nuclides with a relatively large 
Qo value, a correction for the a effect becomes really 
necessary. To reduce the a effect, it is either to develop more 
accurate and precise techniques for a determination or to 
choose the irradiation channels with the a value low enough. 

In the case of the epicadmium neutron activation, Eq, (13) 
can be changed into: 



Z p (a) = |-a(a>(E/)- ai ln(E n | 



(14) 



The values of the error propagation of a in this case were 
carried in Table 8. In this case, it clearly shows the inaccuracy 
of the approximation expression in [6] (Eq. (4) in this report). 
Really, acording to Eq. (4), the error propagation function of 
a in the irradiation with cadmium can be written: 

Z p (a)=|-a(ln(E/)-ln(E7)| (15) 
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TABLE VII 

CALCULATED RESULTS OF Zp(a) FOR CHOOSEN NUCLIDES 



Nuclide 


7-1 channel 


Neutron trap 


Channel 17 




(Dalai reactor) 


(Dalat reactor) 


/ T 1 T T Ti 1 T T C 1 

(1 HE 11a 

reactor) 


45 Sc 


0,03571 


0,01684 


0,02241 




0,01755 


0,01022 




94 Zr 


0,03379 


0,01083 


0,02489 


185 W 


0,02154 


0,00916 


0,01390 


1,7 Au 


0,000 


0,000 


0,000 


98 Mo 


0,11340 


0,07016 


0,07294 


% Zr 


0,19924 


0,14750 


0,12751 




TABLE VIII 




CALCULATED RESULTS OF Zp(a) FOR THE NUCLIDES IN ENAA 


Nuclide 


7-1 channel 


Neutron trap 


Channel 17 




(Dalat reactor) 


(Dalat reactor) 


(THETIS 
reactor) 


45 Sc 


0,00835 


0,00626 


0,00531 


59 Co 


0,11817 


0,08863 


0,07520 


94 Zr 


0,2920 


0,2190 


0,18584 


186 W 


0,05509 


0,04132 


0,03506 


""Au 


0,000 


0,000 


0,000 


98 Mo 


0,1640 


0,1230 


0,10439 


% Zr 


0,17969 


0,13477 


0,11435 



Eq. (15) is different to Eq. (14) by the correctional 
coefficients a ; . However, the value of the error propagation 
function in channel 7-1 of Dalat reactor, for 45 Sc is 0.0083 
from Eq.(14) and 0.2997 from Eq. (15). If the uncertainty of a 
in experiment is 100%, the contribution of uncertanty of a on 
NAA result is 0.83% and 29.97%, respectively. It differs by a 
factor of 30 (!). Similarly, in channel 17 of Thetis reactor, the 
error propagation function for 45 Sc is 0.0053 and 0.1907. The 
difference is huge. This comment is also correct for nuclides 
ha ving Qq<1. It once more confirms that the approximation 
expression in papper [6] is not good for nuclides having Qo<l- 
From Eq. (13) or Eq. (14), we easily estimate the influence of 
a on NAA results, if we know uncertainty of a in the 
irradiation channel. However, for ENAA (epicadmium 
neutron activation analysis) the situation is much more 
dramatic, especially for nuclides with low Q value. 



3) Influence of Uncertainty off on NAA Results 

The error propagation function Z p (f) can 
written: 



be 



z p (f) = 



Q 0i (E I1 )--Q M (E I| r a 



(Qo, +f(E ri r ia )(Q ()I +f(E r ,r 



(16) 



The values of the error propagation function of f in the 
channels; 7-1 and neutron trap of Dalat reactor and channel 17 
of THETIS reactor for the 

nuclides; 45 Sc, 59 Co, 94 Zr, 186 W, 197 Au, 98 Mo, 96 Zr were carried 
in Table 9. The uncertainty of f contributes on the error of 
NAA results is: 



^ = Z(f)^ 

p pW f 



Generally seeing, the uncertainty of f in experiment is 
about less than 4%, therefor, from Table 9, the contribution of 
the uncertainty of f on the error of NAA result is about less 
than 2%. 

TABLE IX 

CALCULATED RESULTS OF Zp(f) FOR THE CHOOSEN NUCLIDES 



Nuclide 


Channel 7-1 


Neutron trap 


Channel 17 




(Dalat reactor) 


(Dalat reactor) 


(THETIS 








reactor) 


45 Sc 


0,512 


0,321 


0,494 


59 Co 


0,398 


0,2697 


0,393 


94 Zr 


0,204 


0,167 


0,225 


185 W 


0,0202 


0,0208 


0,0252 


""Au 


0,0 


0,0 


0,0 


98 Mo 


0,282 


0,323 


0,281 


% Zr 


0,413 


0,566 


0,428 



4) Collective Influence of Uncertainties of a, Eri and f on 
NAA Results 

In view of the above, we can estimate the influence of the 

uncertainties of a, Eri and f on final NAA results. The 
contribution of these parameters on the errors of the analysis 
results is written as: 



Ap 
P 









2 


f 




2 








(Ap 


+ 


Ap 




+ 




a,E r ,a ^ 


I P 


a J 




v P 


e7 J 




v P 


f ) 



(18) 



(17) 



However, as discussion above, the E r effect is negligible 
and can be omitted in Eq. (18). Thus, the contribution on error 
of NAA results in this case is primarily due to the 
uncertainties of a and f. Finally, as well as estimation above, 
this overall contribution of a and f is about 2% on the error of 
NAA results. It was also confirmed by actual analysis. 

IV. CONCLUSION 

For a in the iradiation position relatively small (a «1), 
Eq.(5) is a good approximation to estimate influence of the 
neutron flux characteristics on NAA result using the k - 
standardization method. From this approximative expression, 
the error propagation functions of the parameters were 
presented. They can be used for the estimation of the errors on 
NAA due to the uncertainty of the neutron flux characteristic 
parameters at the irradiation position. From the results of this 
report, it was also confirmed that the approximation in papper 
[6] is only acceptable for the nuclides having Q 0l >1, but not 
for the nuclides having Q 0l <l. 
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